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Abstract 

We discuss systematically several possible inequivalent ways to describe the dynamics and 
the transition probabilities of a quantum system when its hamiltonian is not self-adjoint. 
In order to simplify the treatment, we mainly restrict our analysis to finite dimensional 
Hilbert spaces. In particular, we propose some experiments which could discriminate 
between the various possibilities considered in the paper. An example taken from the 
literature is discussed in detail. 
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1 Introduction 


In ordinary quantum mechanics one of the fundamental axiom of the whole theory is that the 
hamiltonian H of the physical system must be self-adjoint: H = //' . This condition, shared 
also by all the observables of the system, is important since it ensures that the eigenvalues of 
these observables, and of the hamiltonian in particular, are real quantities. However, this is not 
a necessary condition, and in fact several physically motivated examples exist in the literature 
concerning non self-adjoint operators whose spectra consist of only real eigenvalues. 

However, H — /F has an extra bonus, since the time evolution deduced out of H is unitary 
and, being so, preserves the total probability: if T(t) is a solution of the Schrodinger equation 
then ||T(£)|| 2 does not depend on time. This is clear since T(t) = e _ *^'I/(0), 
and since U t = e~ lHt is unitary. Of course, this is false if H F H', and in fact, in this case, 
||T(£)|| 2 does indeed depend on time, in general. Sometimes this is exactly what one looks 
for: in many simple systems in quantum optics, for instance, non self-adjoint hamiltonians are 
used to describe some decay, so that there is no reason for the probability to be preserved in 
time. Other times, one would prefer to avoid any damping, so that the aim is to find some 
way to recover unitarity even when H F iF. This is particularly interesting for people in the 
PT-community, who quite often work with hamiltonian operators which are not self-adjoint, 
but simply pseudo-symmetric or PT-symmetric, IB El. and in fact several attempts have been 
proposed along the years by different authors to discuss this and other aspects of time evolution 
for systems driven by non self-adjoint hamiltonians. Here we refer to [3j-[T^, and references 
therein. However, in our opinion, much more can be said, and using a rather general approach. 
This is exactly what we will do here, in the next section, considering the cases in which the 
eigenvalues of H are all real and commenting on the situation in which some eigenvalues are 
complex. 

In all this paper we will work with finite-dimensional Hilbert spaces. This has two nice 
consequences: the first one is that all the operators involved are bounded (hence, everywhere 
defined) and the inverse, when it exists, is bounded as well. In fact, we are dealing with matrices. 
The second consequence is that we can easily, quite often, discuss examples in terms of pseudo- 
fermions (PFs), TlJ TUI. as we have already recently shown in [13] . We should stress that, 
contrarily to what often stated in the literature, going from a finite to an infinite dimensional 
Hilbert space is an absolutely non trivial task. Therefore, most of our claims, though giving 
indications also in this latter case, are rigorously true only in the present, finite-dimensional, 
settings. We will comment more on this aspect all along the paper. 
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This article is organized as follows: 

In the next section we discuss the general functional structure associated to a non self- 
adjoint hamiltonian, and its dynamics. We also comment briefly on the case of non purely 
real eigenvalues and on finite temperature equilibrium states. In Section 3 we propose different 
definitions of transition probability functions, and we discuss a possible strategy to discriminate 
between them. This is, in fact, the core of our paper since it could be used, in principle, to 
deduce which are the correct Hilbert space, scalar product, norm and adjoint, or, more explic¬ 
itly, which definitions reproduce the experimental data. This proposal is made more precise 
in Section 4, with the aid of an explicit example, originally introduced in [4] and discussed 
here adopting a simple and general pseudo-fermionic representation. Section 5 contains our 
conclusions. To keep the paper self-consistent, we list some definitions and results on PFs in 
the Appendix. 

2 A general settings for H ^ H' 

As we have already said, in this paper we will focus on the easiest situation, i.e. on finite 
dimensional Hilbert spaces. In this way our operators are finite matrices. The main ingredient 
is an operator (i.e. a matrix) H, acting on the vector space C N+1 , with H A H 1 and with 
exactly A^ +1 distinct eigenvalues E n , n = 0,1, 2,..., N. Here, the adjoint W of H is the usual 
one, i.e. the complex conjugate of the transpose of the matrix H. Because of what follows, 
and in order to fix the ideas, it is useful to remind here that the adjoint of an operator X , X', 
is defined in terms of the natural scalar product (.,.) of the Hilbert space V. = (C JV+1 , (.,.)): 
(Xf,g) = (f,X*g), for all f,g E C N+1 , where (f,g) = Ylk=o fk9k, with obvious notation. We 
will consider separately the case in which all the eigenvalues E n are real and the situation in 
which some are complex. In both cases we will assume that each E n has multiplicity one. 

Before starting, it is necessary to clarify some notation adopted in this paper: we will use 
C N+1 any time we want to stress the nature of vector space of our vectors. When it is important 
to stress the topological (i.e. the scalar products and the norms) aspects of this set, we will use 
H instead of CA 44 (and, later, or Hq,). Before starting with our analysis, it is surely worth 
stressing that, with a different language, some of the results discussed in Section 2 can be found 
in the literature, see gusm ns nu for instance. We have decided to include these statements 
here for several reasons: first, they are useful to fix our notation. Secondly, some of the proofs 
discussed here are different, or cannot be found, in the existing literature. Last but not least, 
we want to keep an eye to possible extensions of our results to the situation in which infinite 
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dimensional Hilbert spaces are needed, where unbounded operators most probably appear, with 
all their delicate mathematical aspects. 

2.1 All the eigenvalues are real 

We assume here that H has N +1 distinct real eigenvalues, corresponding to N +1 eigenvectors 
ip k , k = 0,1,2,... ,N: 

Hip k = E k (p k . (2.1) 

The set Jy, = {(p kl k — 0,1, 2,..., N} is a basis for C^ 1 * * * , since the eigenvalues are all different. 
Then an unique biorthogonal basis of Ti, J 7 ^ = {Tfc, A: = 0,1, 2,..., N}, surely exists, [T5l fl6j : 
(<p k ,^i) = S k ,i, for all k, l. Moreover, for all / £ Ti,. we can write / = J2k=o /) = 

J2k=o (^k,f)<Pk- Incidentally, this means that both T v and T ^ are complete (or total): if 
/ G % is such that (tp k , f) = 0, or (T fc , /) = 0, for all k, then / = 0. 

What is interesting for us is that the set J 7 ^ is automatically a set of eigenstates of iU with 
eigenvalues E k : 

H^ k = E k V k , ( 2 . 2 ) 

k = 0,1,2,..., iV. This follows from the completeness of and from the following equality, 
where k is arbitrary but fixed: 

{{H^ k — -Efcdq.) , </?/) = (ifh V k , (pi)—(E k ^ k , <pi) = (T fc , Hcpi)—E k (T fc , (p t ) = (. E t — E k ) 5 kyi = 0, 

for all l — 0,1, 2,..., N. Then (12.2p follows. 

Using the bra-ket notation we can write J2 k=0 = J2k=o = 1, where, for 

all /, g, h e "H, we define (\f)(g\)h := (g,h) f. We introduce the operators S, p = J2k=o \ l Pk)( ( Pk\ 
and Sy = Ylk=o as in |H3- These are bounded positive, self-adjoint, invertible 

operators, one the inverse of the other: Sq, = S'" 1 . We want to stress again that, in our present 
settings, there is absolutely no problem with the domains of these (and other) operators, while 
in H33 we have discussed what happens for infinite-dimensional Hilbert spaces. This passage is 
absolutely non trivial, but will not be considered hereQ. Using standard techniques in functional 
analysis, or direct matrix computations, we can introduce the positive square roots of Sq, and 

1 Operators of this kind have been introduced by several authors, with different names, in recent years. 

Some references are 13 mug, where other references can be found. However, only recently the relevance of 

unbounded operators in this context has been recognized, [18j 19, .2Cjl[2T]. As it is well known, in this case 

several problems, mainly related to domain problems, must be taken into account. 
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S p , and again we have S^ 


1/2 


0 - 1/2 

Ocp 


Useful (and well known) properties are the following: 


S v V n = <p n , Svp n = as well as S 9 H = H j Sq,, S P H ] = HS V . (2.3) 

If we now define H 0 := S\j 2 3 HS p 2 and e k = S^ 2 p k = S p '^ k , k — 0,1, 2,..., N, we see that 

Ho = Hi = S 1 / 2 HKs!J'\ and that £ = {e k , k = 0,1, 2,..., N} is an orthonormal (o.n.) basis of 

% of eigenstates of H 0 : H 0 e k = E k e k . 

Similarly to what is done in many places in the literature, these operators can be used now 
to define new scalar products in C Ar+1 : 


(/, 9>„ := (S/ 2 /. Spg) = (S v f, g), (/, <,}* := (s'Pf, S'fg) = (S*f, g), (2.4) 

for all /, g G C N+1 . Due to the properties of Sq, and S pi these are scalar products, everywhere 


defined on C iV+1 . Of course, the related norms 
instance we can check that 

1 

< 


and 


El 


U' 


are all equivalent!/]. For 


k 1/2 i 


< lis/ 2 ll ll/ll. 


for all / G C N+1 . Then, from a topological point of view, H, T~L P := (c N+1 , (., .)^j and 
Hq, := (C w+1 , (., .) 4 ,) are all equivalent. However, they are different under other aspects, as we 
will show in a moment. In particular, the first difference is in the definition of the adjoint of 
the operators, which is f in T-L, but which becomes b in H P and ft in Hq,: ( Xf,g) = (f,X^g) 
and ( Xf , g)^ = (/, X^g} , for all /, g G CD 44 . It is easy to see that Jf and b are really adjointtjj, 
and to deduce the following relations: 


X* = SvX*S v , X^SpX'Sv and X b = S^S^, (2.5) 

for each operator X on C^ 41 . It is now an easy computation to check that H = H\ and 
that iff = (fU) b . Indeed we have (. Hf,g )* = (f,Hg) 9 and (H*f,g) v = (f,H*g) v , for all 
f,ge C N+1 . 

Remark:— the equalities in (j2.5j) cannot be extended easily if dim('H) = oo. The reason is 
the following: if, for instance, X^ and S ip are unbounded, taken / G D(S,p), the domain of S p , 
there is no reason a priori for S p f to belong to D(X^), so that X^S p f needs not to be defined. 
Therefore, when dim('H) = oo, the three Hilbert spaces are different not only topologically, 
but also as sets. Of course, this cannot happen if dim('H) < oo, since all the operators can be 
defined in all of Ti. 

2 This means that, if a sequence of vectors /„ € C w+1 converges in ||.||, it also converges in ||.||<p and in ||.||^. 

3 For instance, (Xty = X and (XY)$ = Y^X$, for all operators X and Y on C w+1 . 
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2.1.1 The dynamics 

The aspect we are interested in here is the dynamics associated to H and to , since there 
is no reason a-priori for the Schrodinger equation to hold also in this case. This problem has 
been discussed by many authors over the past decade, see 13 El El HD] for instance, with a 
particular attention, sometimes to the unitarity of the time evolution. However, our point of 
view is slightly different since we are more interested in: (1) motivating the extended Schrodinger 
equation usually taken for granted in the literature, and (2) trying to understand if it is possible 
to deduce which is the best choice (with best to be somehow defined!) for the time evolution of 
a system, in this case. 

Our starting assumption is that, since H 0 is self-adjoint, it produces a standard Schrodinger 
equation for the wave function $(t): iQ(t) = H 0 Q(t), whose solution is $(t) = e~ lHot Q(0). 
This function satisfies the following equality: p{t) := ||$(t)|| 2 = ($(t),$(t)) = p(0). Then, 
probability is preserved. What is interesting is that, using our results, and in particular the 
similarity equations H 0 := S\j 2 HS^J 2 = S^J 2 S ^ 2 , two related Schrodinger-like equations 
naturally arise, with simple computations, from the one for <3>(t): 

i4>^(t) = H$^(t), and i&^t) = H^ v {t), (2.6) 

where <f>^(t) = S'^ 1 ' /2< h(t) and $ v (t) = S'^ 1 ^ 2< f>(t). Of course, both these equations can be easily 
solved: 

$*(t) = e" im <M0) = e- iHt S~ 1/2 <S>( 0), %(t) = e - iHU $ v (0) = e~ iHU 0). (2.7) 

Now, the equations in (12.6(1 suggest that H and produce, by themselves, two Schrodinger- 
like equations which could be, therefore, taken as the starting points to describe the dynamics 
of the given system. This is, in fact, what it is done, quite often, in the literature, 13 ei m 
etc.: one uses a non self-adjoint hamiltonian H , and simply writes the equation of the dynamics 
for the wave function £(i) as in (12.6(1 : i£(t) = H£(t). The two equations in (12.6(1 . which are 
deduced directly from i$(i) = H 0 $(t), suggest that this is in fact reasonable. Of course, the 
old probability is not preserved for <f>,j,(t) and ^(t), at least if we stick with the norm ||.||. 
But there is no a-priori reason to adopt this particular choice. One could use 11.11^ or ||.||^, and 
indeed we have 

Pv(t) := ||<Mt)|| 2 = (<Mt), <M t)) v = fv(0) 

and 

Pq/(t) := ||$^(t)||^ = ($^r(t), $^r(t))^ = P'I'(O), 
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for all t e 1. Then we see that, if for some reason we are interested in having the probability 
associated to preserved in time, the natural framework to use is not the one provided by 

H but that given by H ip . We will go back to this point in Section 3. 

This freedom of choice is reflected by the dynamics of the operators in the Heisenberg 
representation: in standard quantum mechanics, i.e. working with H 0 = Hq in H, the recipe 
to deduce the time evolution of the observable X is simple: its time evolution, X(t), should 
have on the vector £ describing the physical system at t = 0 the same mean value of X on the 
time-evoluted vector £(t) = e~ lHot ^\ 

(fWWf) = 

Then X(t) = e !fl °We _!floi . If we repeat the same steps, replacing Ho with if, we deduce that 

X n (t) = e iHh Xe~ iHt . (2.8) 

But this is indeed not the only possibility, and it is not even the most convenient, since during 
this time evolution we lose three crucial characteristics of the standard time evolution: (1) the 
operators e lHtt and e~ lHt are no longer unitary. Sometimes, however, this is exactly what one 
looks for. As already stressed, this is what is done, for instance, in quantum optics, to describe 
some damping, [22] [23j; (2) An annoying consequence of (12.81) is that it is not so easy to find 
integrals of motion for the system, since [H,X] = 0 does not imply that X(t) = X(0). Then, 
it is no longer true that the observables which commute with H do not evolve in time; (3) A 
more serious difficulty is that the time evolution is no longer an automorphism of the set of 
observables, since in general (XY)%{t) ^ Xy_(t)Yy_(t), and this complicates in an enormous 
way all the computations. For instance, in general, it is absolutely non trivial to deduce the 
analogous of the Heisenberg equation of motion for X(t). 

However, there is a possible way out, and it has again to do with a suitable choice of the 
scalar product. In fact, if we work with Hq, rather than with "H, which is a natural choice since 
H is self-adjoint in Hy, we can define a different map: t —> X^(t) working as before: 

(Zv(t),X£is,(t)^^ = <^£*, A«*(t)£tf^ , 
and the result is the following: 

X H „(t)=e iHt Xe- iHt , (2.9) 

and problems (1), (2) and (3) are solved, paying the only price to work in Hy, i.e. to replace 
(.,.) with (.,.)^. Apparently, this is not a big price, indeed. This is not yet the end of the 
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story, since one may wonder why H should be better than ifh In fact, there is essentially no 
difference, a priori. In fact, if for some reason we need to use W rather than H, the choice of 
working in "H would create the same kind of problems (1), (2) and (3) as before. The possible 
way out is now clearly to use as the natural Hilbert space to work with. In this case, the 
time evolution of X looks like 

X Hv (t)=e iHh Xe~ iHh , (2.10) 

and again those problems are solved. 

Summarizing, we can say that, under our assumptions on the E n , working with a non self- 
adjoint hamiltonian H gives a lot of freedom: we have three Hilbert spaces, with their scalar 
products, their norms and their involutions, and we further have several possible definitions 
of the time evolution of any observable, each one with pros et contra. Of course, we are left 
with a very natural, and deep, question: how should we choose the right Hilbert space? Or, 
even better, who or what decides what is right ? We will suggest a possible way to answer 
these questions in Sections 3 and 4. This is a problem which, in our knowledge, has not been 
considered in details in the literature so far. 

Remarks:— (1) Quite interestingly, and as it is implicit in the work of several authors, 
see ME] among others, in this context it is the hamiltonian H itself which somehow fixes 
its preferred Hilbert space! This is because both (.,.) and (., .)^ are defined via S v and S 
which are deduced, in turns, by the eigenvectors of H and Hh This is similar to what happens 
in algebraic quantum dynamics, see |24| and references therein, where the hamiltonian (self- 
adjoint, in that context) is used to define a suitable topology on the algebra of the operators 
needed in the description of the physical system. 

(2) Of course, the possibility of having different definitions for the time evolution is reflected 
in the definition of the equilibrium states for non-zero temperature. More explicitly, if we use 
Hq to define the dynamics as X(t ) = e lHot Xe~ lHaT , the natural choice of equilibrium state is 
the Gibbs state 

wo(X) := T( r (e-^X), 

where Zq = tr(e~^ H °) /3 = ^t, T being the temperature of the system, and tr(A ) is the trace of 
the operator A. Here we are fixing to one the Boltzmann’s constant K. Of course, this is not 
the most reasonable choice if one imagine that the time evolution of the observable X is given 
by X H<j! (t) = e lHt Xe~ lHt or by X Hv (t) = e lHit Xe~ iHtt . To these choices, in fact, it is more 


natural to associate the following states 

u> 9 (X) := 2-tr{e-e H X), u> v {X) := ±-tr(e~^X), 

Zvp Zi^p 

where = tr(e~^ H ) and Z v = tr(e~^ H ). To avoid mathematical difficulties, we are thinking 
as before that the Hilbert space of our system is finite dimensional, so that all the quantities 
introduced here are well defined. It is not hard to check that all these states satisfy a KMS-like 
equilibrium condition like, for instance, 

= ^(BA^t + f/3)), 

for all observables A and B and for all possible choices of j}, i.e. whenever we use the same 
hamiltonian to define both the time evolution of the system and the state. Hence, loosing 
self-adjointness does not imply many changes from this point of view. However, it could be 
interesting to check what happens adopting, for instance, H 0 to define the time evolution and, 
a different but somehow related operator H , to define the state. 

2.2 Not all the eigenvalues are real 

Let us now briefly consider what happens when we abandon the assumption that E k is real for 
all k. This might be relevant for including, in our scheme, some effective hamiltonians used 
in different, usually non conservative, contexts, see for instance |22[ |23lj for an application to 
quantum optics. In this case, equation Q 272]) must be replaced by 

= ( 2 . 11 ) 

k = 0,1, 2,..., N, which might seem to be a minor difference. This is not so. In fact, in this 
new situation, H and are no longer isospectral, and for this reason no invertible intertwining 
operator exists between H and Hf This makes the framework of the system a little poorer 
than before. In fact, the existence of intertwining operators proved to be quite important in 
many physical systems, as discussed in many papers on this subject, Even more: suppose 
H f and suppose that 5s(E n ) f 0 for at least one n. Then it is easy to show that, contrarily 
to what happens for real eigenvalues, there exists no scalar product -C .,. on C Ar+1 such 
that <C Hf,g 3>=-C f,Hg ^>. The proof is trivial, but in our opinion is worth giving: since 
•C H(p n ,ip n ^>= E n -C (pniVn and (p n ,Hip n »= E n ip n ,(Pn >•, it is clear that, if 
E n E n , Hip n ,ip n tp n ,Htp n 3>. Hence H is not self-adjoint with respect to this 

different scalar product. The conclusion follows from the arbitrariness of <C .,. 3>. 
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Then what we did before is no longer true in this new settings: more in details, Jy, and J A 
are still biorthogonal bases (or, most probably, biorthogonal sets if dim('H) = oo), and S v and 
Sy are again one the inverse of the other, positive and self-adjoint. Their square roots exist, 
and e n = Sjj/Vn is still an eigenvector of H 0 = S',]/ 2 HS l J 2 with eigenvalue E n . Also, £ = {e n } 
is a basis for E. However, H 0 A Hj } , S^H A and Hle n = E n e n . In particular, this 

shows that Sq, is no longer an intertwining operators between H and h/A 

A consequence of these facts is that, assuming as before = if$^(f) and i<&<p(t) = 

H^,p(t), their solutions do not preserve probabilities. So, in this case, there is no leading rule 
to follow, apparently. 

Remark:— It is clear that when E is finite dimensional, the presence of complex-valued 
eigenvalues makes not a real big (technical) difference: in fact, we are simply dealing with 
(N + 1) x (N + 1) matrices! However, apart from what we have already discussed, we can 
expect more problems when E is infinite dimensional. The analysis of this situation is work in 
progress. 


3 Transition probabilities and consequences 

In this section we propose three different definitions of transition probabilities, and we deduce 
the different results which are obtained out of these definitions, suggesting some experiments 
which make it possible to discriminate among them. In this way we should be able to de¬ 
cide which is the correct framework to adopt when non self-adjoint hamiltonians, with real 
eigenvalues, are those who naturally describe a physical system S. 

We first briefly recall what happens in the standard situation, i.e. when the dynamics of 
S is driven by a self-adjoint hamiltonian Ho, with an o.n. basis T e = {e k , k = 0,1,..., N} of 
eigenvectors: H 0 e k = E k e k , Vfc. In this case, E k e M automatically, for all k. If $0 = Ylk =0 Cfcefc 
is the state of the system S at t = 0, c k = (efc,<ho), then its time evolution is clearly given 
by <3 ){t) = e~ lHot ^o — Yhk=v c k&~ lEkt &k- Here it is not a major request assuming that <R(t) is 
normalized for any t, since e~ lHot is an unitary operator. In order to compute the transition 
probability to a final state, described by the normalized vector $/, we just need to compute 

(3.1) 

which, because of the Schwarz inequality, is clearly always between zero and one: G 

[0,1], for all t. Of course, this is strongly related to the fact that both $ f and <3>(f) are normalized 
in E, and this is possible for all teR since Ho is self-adjoint. Just as an introductory example, 
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let us now see what happens if $o = <^o and $/ = 'ho, where ipo and 'ho are two vectors of 
two generic biorthogonal sets = {p n } and = {'h n } ) none of which made of eigenstates 
of H 0 . Of course, we are interested in this particular situation in view of our next extension 
to non self-adjoint hamiltonians. In this case, since normalization is decided by the condition 
( 990 ,'ho) = 1, there is no guarantee that ||</?o|| = ||'ho|| = 1. For this reason, formula (13.11) 
produces 




11^0 


I'hr 


|('h 0 ,e-^Vo>r = 


ll^o 


I 'hn 


N 


^ d k p k e 


r), 0 iEkt 


k =0 


(3.2) 


where d k = {e k ,Po) and p k = (efc,^/ 0 ). This formula shows that, in particular, P$ 0 ^<s> f (t) does 
not depend on time if all the pairs (d k ,p k ), except one, are zero. Also, if the eigenvalues E k 
are commensurable, is a periodic function. Not surprisingly, if tp 0 = 'h 0 , we get 

(d) = 1- 


Let us now consider the non self-adjoint settings described in Section 2. Hence we have 
H 7 ^ H\ and two sets = {p n } and = {'h ra } of biorthogonal eigenvectors of H and 
ifh Hp n = E n ip n , if t 'F n = E n ty n , {(p k , T n ) = 5 ktU . Notice that here we are restricting to real 
eigenvalues, E n G M, as we will do in all this section. 

As we have discussed before, it is natural to assume that the wave function of the system 
S described by H satisfies the equation *$(£) = with initial condition <b(0) = <h 0 . 

Then, since the <pn s are eigenstates of H , the natural choice to compute <L(f) is to expand this 
unknown function in terms of Hence we get 

N 

$(t) = e ~ lHt ^ 0 = yc k e~ lEkt (p k , (3.3) 

k =0 

where c k = ('h*,, $ 0 ). It is clear that e~ lHt is no longer a unitary operator, so that there is no 
reason for || < I ) (t)|| to be always equal to one, even when j| < I ) o|| = 1. For this reason, at least if we 
use the norm ||.||, a time dependent normalization must necessarily appear into the game. As a 
matter of fact, more than one such normalization will now be considered, each one related to a 
possible different definition of the transition probability in the present context. We should also 
stress that the necessity of introducing some normalization in similar contexts was already clear 
to many authors, already several years ago, [26], and recently reconsidered by other authors, 
[3], as well as by many others. 

In practice, going from self-adjoint to non self-adjoint hamiltonians opens the possibility of 
having several possible (apparently) inequivalent definitions of transition probabilities, all of 
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which appear to be absolutely reasonable. The ones we introduce here are the following: 

<*/,*(*)>„ 2 

*/IUI*(*)L ’ 

(3.4) 

which are the most plausible definitions, but not the only ones, as we will discuss later. They 
look different since they involve different scalar products and different norms. We will see that 
there is more than this: they really produce different results, so that they are not physically 
equivalent at all and it should be possible, in principle, to discriminate among them, in order 
to understand which is the most appropriate expression of the transition probability, and why. 
Incidentally we observe that, with these definitions, 6 [0,1] for all t > 0. Here 

stands for P£ 0 ^ f or P£ 0 _> $/ . 

Remark:— At a first sight, P$ 0 _t$ may appear as the more natural choice, since H is 
self-adjoint with respect to (., .) 9 , and therefore e~ lHt is unitary with respect to this scalar 
product. However, here we are not considering unitarity of the time evolution as our main 
requirement. We are much more interested in a comparison between our theoretical results 
and some experimental data, and a good agreement is not necessarily ensured by the unitarity 
of e~ %Ht . Notice also that, if we replace H with H\ this unitarity request would suggest, of 
course, to use P§ 0 _>.$ rather than Pjf ^, Pf . But, again, this is not really our main criterion. 

We begin our analysis with a simple situation: let us assume that S is prepared in the 
following linear combination: $o = Va + <fb, where 0<a<N,0<b<N, and a ^ b. Then 
$(t) = e~ lHt $o = e~ lEat (p a + e~ lEbt ifb- Hence we get 

\\®(t )\\ 2 = \\Va \\ 2 + ||rf+ (e l{Ea ~ Eb)t {(p a ,<Pb) +c.c .) , 

ll*(*)ll! = 2, 

and 

= \\<Pa\\l+ \\Ml+ {e l{Ea ~ Eb)t (<p a ,<p b ) v + c.c^ , 

where c.c. stands for complex conjugate. 

If we now take $/ = <fj, for some j = 0,1, 2,..., N, we conclude first that 

= 2 + ^j,b) • ( 3 - 5 ) 

Hence, according to this rule, a transition between <h 0 = <p a + Lp b and <pj is possible if and only 
if j — a or j = b. In all other cases, a similar transition would not be allowed. However, this is 


P 


$0 —f 


(t) ■= 


<$/,$(*)> 




p'V 


(t) ■ = 


(*/,*(*)>< 
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not the same conclusion we get considering the other possible definitions. In fact we find 


Pc 


^O - ^ / 


(*) 


{<Pj, <Pa) \ 2 + I {<Pj, Vb) I 2 + {e l{Ea Eb)t (y?q, <Pj) {(Pj, Lpb) + c.c) 
Ill'll 2 (ll^all 2 + ll^fell 2 + (e i ( Ea ~ Eb '^ t (ipa, (fib) + C.c)) 


and 




in, v*) v I 2 + 1 in, ft) v l J + (« i<£ “ Eb)t {fa, fi) v (fj,fb) v + c.< 
WflWl (il<A>IIJ + \\fb\\l + (fa,fb)„ + c.c 


(3.6) 


(3.7) 


It is clear that, if we perform an experiment on S, and we find that the transition probability 
from a vector <h 0 = <p a + ip b to with j ^ a,b , is different from zero, P$ 0 ^$ / (t) must be 
discharged, while the other two possible definitions could still be correct, in principle. 

In a similar way, if we fix now d> j- — du leaving unchanged d> 0 , we deduce that 


P$o —(.") — 


5 j} a + 8. 


j,b 


Id' ■ 


(Il'T’oll 2 + IMI 2 + (e* (E “ Eb)t {fa, fb) + C.c)) ’ 


(3.8) 


p v I' 


(t) 


and 


<*i, *a} | 2 + | (%, %) | 2 + (e^ Ea ~ E ^ (T a , T,) (^-, %) + c.c) 

2||^||| 


(3.9) 


/ 


(0 


1 {<Pj,<Pa) 

| 2 + 

(<Pj, (fb) 

2 + ( e i{Ea Eb)t {<Pa, <Pj) {<Pj, <Pb) + C.c) 

ll^all 2 

+ 

+ 1 

f e i(E a -E b )t <p 6 )^ + c.c) 

1 


(3.10) 


We see that, with these choices of initial and final states, we should get a zero transition proba¬ 
bility if j ^ a, b, at least if the correct transition law is given by P$ 0 - i .$ f (t). Suppose then that, 
making an explicit experiment and choosing j ^ a , b, we observe a zero transition probability 
from $ 0 to d>j. Then we are forced to assume that the only possible correct expression for 
such a probability is exactly the one given by P<S> 0 ^f (t). Summarizing, the probability func¬ 
tions proposed here all make perfect sense, and only some experiment can discriminate between 
them. 

We will say more in the next section, where an example is discussed in many details. 

Remarks:— (1) It remains open the case in which some of the eigenvalues of H are complex. 
In this case we have already seen that problems may arise in the general settings proposed here, 
so a deeper analysis is required. Notice also that none of the formulas above describe damping: 
only oscillations are allowed! The reason is clear: our analysis, here, is restricted to real 
eigenvalues. In other words, even if the hamiltonian is non self-adjoint, its eigenvalues are real 
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numbers anyhow. In order to get some decay, we need also some eigenvalues with non zero 
imaginary parts. 

(2) In principle, we could repeat the same analysis using W instead of H , but we will not do 
it here. In fact, we do not expect any essential difference. On the contrary, we expect that the 
conclusion would be quite similar: we can still consider three possibilities as in (13.41) . replacing 
H with H\ but only experimental results could say which one is correct. 

(3) Once we have been able to decide which one is the probability transition which is in 
agreement with the experiments, it is clear that we have a preferred scalar product , and our 
suggestion is that this is the one to be used in the computation, say, of the norm of the 
wave-function, in the determination of the adjoint of the operators (and of the observables in 
particular), and so on. 


4 A detailed example 


The system we will discuss here was considered recently in [3], and it is described, in its 
simplified version, by the non self-adjoint hamiltonian 


Hsds — ~9 


0 1 - k 

1 + k 0 


) 


where g e M and k e] — 1,1[, and we have put h = 1 to simplify the notation. Of course, the 
interesting situation is when k ^ 0, since otherwise Hsds — Hsds- This can be written in 
terms of pseudo-fermionic operators, see Appendix, by introducing 




1 1/a 


, and, taking p = —gV 1 — k 2 and to = 2p, 


we 


where a = Then N = ba = \ ( ^ 

deduce that Hsds = uN + pi. The eigenvectors of Hsds, and of its adjoint, can now easily 
deduced: 


Po = 


1 

—a 


ipi = bip 0 = 



To — N, j/ 


1 

— 1/a 


Ti = a t T 0 = 


1 

1/a 


5 


where N v Nq, = to guarantee that (<pk,^i) = 5k,i, k,l = 0,1. Then Hip>k = EkPk, with 
Eq — p — —Ei, which are both real for the range of k allowed. Analogously, we can explicitly 


14 




check that H^k = E^k, k = 0,1. A simple computation shows that El=o \ L P k )( = 

ELo l^fe)(^fe| = A while 


S<p = ^ Vk){Vk\ = 2|AE 


fc =0 


1 0 
0 a 2 


, S* = ^|tt fc )(*fcl = 2|W*i 


fc =0 


0 1/a 2 


Then, as expected, S v = S^, 1 = S and they are positive operators, with obvious pos¬ 
itive square roots S^ 2 and S^ 2 . Due to the fact that these are diagonal matrices, it is 
a particularly simple exercise to check that all the properties listed in the Appendix are 
indeed satisfied: S^k = <Pk, S^N = N^Sy, and so on. Also, we could introduce c = 


q !/2 q 1/2 _ 1 [ 1 !/« 

O^j LLQ(p — 2 


, N 0 = c f c = | 


V2i\yiV*| 


-1 -1 

, and e\ = S^Vi = y/^NplNy 


1 1/a 
1 1 


1 /2 

, and the o.n. vectors eo = (po = 
, as well as a self adjoint hamiltonian 


Hr, similar to H and to h/A 


H 0 = Sl /2 HS 1 J 2 = p(l-2N 0 ) = p 


0 -1 

-1 0 


The different scalar products we can introduce in C 2 , given / = 


/o 

h 


and g = 


9 o 
,9i 


are (/,s) = /o9o + fi9i, (f,g)i t = 2|AA| 2 (f Q g 0 + jpfigi) and (f,g) v = 2|A/^| 2 (/ 0 ^o + « 2 /i9i)- 
Moreover, the new adjoints for H , other than fA, using (I2.5jl are found to be 

(i+D 2 


H b = = 


0 

-9 \ (l-k) 2 

1 +k 


1 -k 
0 


while 


if# = S^fAS* = - 


-9 


0 1 - k 

l + k 0 

Now, it is clear that if# = H. Also, it is a simple exercise to show that iA = (iA)#, as we have 
seen in Section 2 for general reasons. 


Remark:— the use of PFs here could be thought as not really essential. And in fact, it is 
not. ffowever, as we have discussed in |T3], it provides a sort of elegant and unifying language 
for many finite-dimensional systems previously introduced in the literature by several authors. 
All these systems share a somehow common structure: biorthogonal sets, lowering, raising and 
number-like operators, nice anti-commutation rules, and so on, and all these features are quite 
naturally described in terms of PFs. 
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4.1 The dynamics and transition probabilities 


Let us now consider what happens if we assume that H drives the dynamical behavior of the 
system S we are considering in this section, via the equation z<j>(i) = HQ(t). In particular, 
following what we have discussed in Section [3] we compute the different transition probabilities 
for different choices of $o and $/. To avoid useless complications, from now on we fix N, p = 

N * = TI- 

To begin with, let us take <f>o = +o and = Iq. Then <3>(f) = e and we get 
P$ 0 ^<f> s (t) = 0) P®o^<t> f (0 = P®, = p _|_ ^ • 


As we see, these are all constant in time. It is clear that, if in an experiment, we compute the 
transition probability from Lp 0 to 4>i, we could conclude that is the correct definition 

of the probability only if the result of the experiment gives zero. Otherwise, if we get pcpp we 
are not in a position to choose between and since they coincide. Then 

we need a second experiment. In particular, we can repeat the same measure, assuming again 
that $0 = (fioi but asking what changes if we now take $/ = tpi. This is a better choice since 
we get three different results: 


— k , P$ 0 ->$f(t) ~ 0) 




Ak 2 

(P + l) 2 ’ 


so that, if k T 1> a single measure would be enough to discriminate between the three definitions. 

Remarks:— (1) If we take k — 0 all the functions above reduce to zero. This is expected 
since, in this case, H becomes self-adjoint, and therefore <pk coincides with and becomes 
an orthonormal set. Hence, in both cases considered above, we are asking which is the possibility 
that a system, originally prepared in an eigenstate of H , evolves toward a different eigenstate 
of H. Of course, in absence of interactions this possibility is zero, and this is exactly our result. 

(2) The explicit results above show that each transition probability assumes values in [0,1], 
as it should. 


Because of the special conditions considered here, we have deduced probabilities which are 
constant in time. In order to get time-depending probabilities, we need to consider different 
initial conditions on S. Let us consider now $0 = To + +i- Hence $(£) = + e lpt( Pi- Let 

further the final state be $/ = To- Then, with simple computations, we get: 

D 1-fc2 d* u\ _ 1 + k2 + 2k cos(2pt) 

*° ) 2(1 — k cos(2pt)) ’ j 2(1 + k 2 ) 
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and 


The first obvious remark is that, when k = 0, these three functions collapse to the same value, 
The reason is clear, see Remark (1) above. When k ^ 0, the probabilities are indeed different 
and again it should be possible to discriminate between them with some experiment. 

Summarizing, to deduce some internal coherence of the whole framework, we should be able 
to decide first which one, between the different transition probabilities, is the one in agreement 
with the experiments, since they all have essentially the same mathematical properties. This, 
we believe, is crucial for a deeper understanding of the theory. 

5 Conclusions 

After a general discussion on the dynamical problem generated by a non self-adjoint hamilto- 
nian, we have seen some consequences of our choices in the computations of several, inequivalent, 
transition probabilities. With the help of a simple example, we have proposed a way to choose 
the correct settings to be used. This could be useful to shed some light on this kind of systems, 
and in fact this is the core of the paper: as we have already pointed out several times, many 
authors have discussed in recent years the dynamical problem associated to non self-adjoint 
hamiltonians. However, in our knowledge, not many attempts have been made to clarify which 
one, among all the possible definitions, is really compatible with experiments. 

It is clear that the extension of our analysis to infinite dimensional Hilbert spaces is highly 
non trivial. In fact quite often the operators involved turn out to be unbounded. This has 
consequences on the domains of the various operators, on the nature of the biorthogonal sets 
of eigenvectors of H and H ', on the inequivalent topologies which can be introduced in the 
Hilbert space of the system, and so on, im However, the fact that H has a purely point 
spectrum needs not to be true anymore, and this complicates quite a bit the treatment in 
concrete situations. Even from this point of view, the use of PFs appears to be a good choice, 
since it allows a natural extension to infinite dimensional spaces in terms of pseudo-bosons, 
which share with PFs many properties, m 

We should also mention that the case in which some eigenvalues of H have non zero imagi¬ 
nary parts is surely more complicated, especially in view of its extension to infinite dimensional 
Hilbert spaces. The reason has been discussed in Section 12.21 Hopefully, the analysis of this 
situation will be undertaken soon. 
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Appendix: Something on PFs 

The starting point is a modification of the CAR {c, <A} = cc* + c* c = 1, {c, c} = {cR c^} = 0, 
between two operators, c and A acting on a two-dimensional Hilbert space "H. The CAR are 
replaced here by the following rules: 

{a, b} = 1, {a, a} = 0, {b, b} = 0, (1.1) 

where the interesting situation is when b ^ a'. These rules automatically imply that a non zero 
vector, ip 0 , exists in V. such that aip Q = 0, and that a second non zero vector, T 0 , also exists in 
V, such that T 0 = 0, m- In general (^ 0 ^ T 0 . 

Let us now introduce the non zero vectors ipi = hip® and Ti = <RTo, as well as the non 
self-adjoint operators N = ba and W = . We also introduce the self-adjoint operators SR 

and SR via their action on a generic / e V.: 

i i 

Sipf = Syf = J^(T n ,/) T n . 

n=0 n= 0 

Hence we have: 


1 . 

atp i = <p o, & f T i = T 0 . 

2 . 

Ntp n = rup n , 91T n = nT n , 

for n — 0,1. 

3. If the normalization of ipo and To are chosen in such a way that (ipo, To) = 1, then 

R Pki T n ) <5fc,nj 


for k,n — 0,1. 
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4. Sep and Sy are bounded, strictly positive, self-adjoint, and invertible. They satisfy 

im<m 2 +im 2 , ii^ii<ii^ir+ii^ir, 

Sep^&n Pni S\frip n T n , 

for n — 0,1, as well as S v = S jjj 1 . Moreover, the following intertwining relations 

SvN = VIS*, s v m = N Sip ^ 

are satisfied. 

The above formulas show that (i) N and 04 behave essentially as fermionic number operators, 
having eigenvalues 0 and 1 ; (ii) their related eigenvectors are respectively the vectors of Jy = 
{</? 0 ,</y} and Tq, = {T 0 , \&i}; (iii) a and are lowering operators for Jy and T& respectively; 
(iv) b and a t are rising operators for Jy and Tq, respectively; (v) the two sets Jy and Tq, are 
biorthonormal; (vi) the operators S v and Sq, are self-adjoint, bounded, invertible, with bounded 
inverse, and map Jy in Tq, and viceversa; (vii) S v and Sq, intertwine between operators which 
are not self-adjoint. Moreover, see Emua, Tp and Tq, are automatically Riesz bases for %. 

Remark:— For completeness we have to mention the paper by Bender and Klevansky, 
[27], where similar generalized anti-commutation rules were introduced, but with a different 
perspective. 
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